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Abstract. : In this paper, we study complex Wishart processes or the so-called 
Laguerre processes {Xt)t>o- We are interested in the behaviour of the eigenvalue 
process, we derive some useful stochastic differential equations and compute both the 
infinitesimal generator and the semi-group. We also give absolute-continuity relations 
between different indices. Finally, we compute the density function of the so-called 
generalized Hartman- Wats on law as well as the law of To := inf{t, det(Xt) = 0} when 
the size of the matrix is 2. 



1. Introduction 

The Real Wishart process is a symmetric matrix-valued process which was introduced 
by M.F.Bru (1989) as follows : Let Bt = {Bij(t))ij be a n x m Brownian matrix and 
define Xt = BjBt. The process {Xt)t>o satisfies the following stochastic differential 
equation (SDE) 

dXt = BjdBt + dBjBt + nimdt 



XtdNt + dNjy/Xt + nl^dt, Xo = bTBq 



where Im denotes the unit matrix, the superscript stands for the transpose, y/Xl is the 
matrix square root of the positive matrix Xt and {Nt)t>o is a m x m Brownian matrix. 
This process is called the Wishart process of dimension n, of size m, starting from Xq 
and is denoted M^(n,m,Xo). Then, the W{8,m,XQ) where 5 runs over the Gindikin 
ensemble {l,...,m — l}[j\m — l,oo[ is defined as the unique solution of the latter SDE 
with 6 instead of n. Thus, it can be viewed as an extension of the squared Bessel process 
to higher dimension. In this way, Donati et al. (2004) tried to extend some well known 
properties of the squared Bessel processes to the matrix case and derived expressions 
such as the Laplace transform and the tail distribution of some random variables, in 
which many multivariate special functions of symmetric matrix argument appear, such as 
Gamma, modified Bessel and hypergeometric functions (Muirhead 1982). However, the 
latter is quite complicated to deal with and to our knowledge, there are no more precise 
results on the law of these variables. Nevertheless, in the complex case, hypergeometric 
functions of Hermitian matrix argument can be expressed as a determinant of a matrix 
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whose entries are one-dimensional hypergeometric functions. In fact, Gross and Richards 
(1991) estabhshed the following result : 

det{x'l~\^q{ai-j + l,...,ap-j + l,...,bq-j + l;xi)) 
pPq{ai, . . . ,ap,Oi, . . . ,Oq;Ji. ) — V{X) 

where X is a m x m Hcrmitian matrix, (xj) are its eigenvalues, denotes the stan- 
dard hypergeometric functions with scalar argument, V{X) = Yl^^jixi — xj) is the Van- 
dcrmonde determinant and pFq is the hypergeometric function with Hermitian matrix 
argument defined by : 

pFq{ar,...,apM,..., bq, ^) " 1. (^^)^ . . . (5^)^ k\ 

where r = (ri, . . . ,Tm) is a partition of length < m and of weight k (i.e. ti > T2 > 
■ ■ ■ > Tm,J2i^i = ('*)t is the generalized Pochammer symbol and Cr is the so-called 
zonal polynomial. We refer to Macdonald (1995) for further details and Lassalle (1991), 
(1991a) for analogous expressions for multivariate orthogonal polynomials. The dctcr- 
minantal representation above is due to the fact that the zonal polynomial is identified 
with the (normalized) Schur functions defined by : 

, , det(xp+"--^) 

Sr[X,,...,Xm)- 

Consequently, one can use integral representations as well as other properties of stan- 
dard hypergeometric functions to get, at least when m = 2, some results that are till now 
unknown in the Wishart case. The rest of this paper consists of seven sections, which 
are respectively devoted to the following topics: in section 2, we introduce the Laguerre 
process of integer dimension. In section 3, we study the behaviour of the eigenvalue 
process. Then, in section 4, we define the Laguerre process of positive real dimension. 
Section 5 is devoted to the absolute-continuity relations, from which we deduce the 
Laplace transform of the so-called generalized Hartman- Watson law as well as the tail 
distribution of Tq, the first hitting time of 0. In section 6, we focus on the case m = 2 
for which we invert this Laplace transform, and finally, in section 7, we compute the 
density of 5*0 := l/(2ro). 

2. Laguerre Process of integer index 

Let B be a nx m complex Brownian matrix starting from Bq , i.e. B = (Bij) where the 
entries Bij are independent complex Brownian motions, so we can write B = B^ + iB^ 
where B^,B^ are two independent real Brownian matrices. We are interested in the 
matrix- valued process Xt := B^Bt which satisfies the following SDE : 

dXt = dB^Bt + B^dBt + 2nlmdt (1) 

Definition 1. {Xt)t>o is called the Laguerre process of size m, of dimension n and 
starting from Xq = BqBq, and will be denoted by L{n, m,Xo). 

Remark 1. For m = 1, {Xt)t>o is a squared Bessel process of dimension 2n, denoted by 
BESQ{2n,Xo). 
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Remark 2. Set Xt = {Xij{t))i^j. One can easily check that 



dXii{t) = 2y/Xii{t)d-fi{t) + 2ndt 1 <i <m, 
where (7j)i<j<m are independent Brownian motions, thus, Xa is a BESQ{2n, (Xo)i 
Remark 3. The equation above imphes that : 



d{tv{Xt)) = 2y/t^{Xt)d(3t + 2nmdt (2) 

where /? is a Brownian motion. Consequently, (ti(Xt))t>o is a BESQ{2nm, tr{XQ)) 
of dimension 2nm starting from tr(Xo). One can also deduce from ([1]) that for every 
i,j,k,l G {!,••• ,m}: 

{dXij,dXki)t = 2{Xii6kj + Xkj6ii)dt, 

which differs from equation (1-1-5) derived by Bru (1989) since for a complex Brownian 
motion 7, we have ^(7,7)^ = and (7,7)4 = 2t . 

2.1. Infinitesimal generator. Let Hm,Hm be respectively the space of m x m Her- 
mitian matrices and the space of m x m positive definite Hermitian matrices. On the 
space of Hermitian matrix-argument functions, we define the matrix-valued differential 
operators : 



d /' d \ d / d \ d f d d 

i- 



dx ydxjkjjf^ dy \dyjkj dz \dxjk dyjk 

We also define: 

dV ^ d'^ f d d\ ^ 52 



dz 



■.= y^ , — — =y 

J ij dzikdzkj' \dxdyj-j ^ dxi^dykj 

Proposition 1. Let functions f satisfying: 

df df df df r „ ■ ■ 

a — = ^ — ' a — = ~^ — for all 
dxij dxji dyij dyji 

Then, the infinitesimal generator of a Laguerre process L{n,m,x) is given by: 

d\. . , ^ f dV^ , ^ f d^^ 



££ = 2n tr(5R J ) + 2[tr(x3f? J ) + tr(y9 J )] (3) 
d 

where — is the operator defined above. 
dz 

Remark 4. Using the fact that xJ- = x, y'^ = —y and iT{AB) = tr(i?^) = iT{B'^ A^) for 
any two matrices A and B, we can see that 

tr ( ——\ - tr (—— \ -tr( ——\ ^ ( ^{ — ?\ - 2tr ( — — 
^ \ dy dx J ^ \dx dy^ J ^ \ dx dy ) \^ dz ) ^ \ dx dy 



3. Eigenvalues of Laguerre Process 



In this section, wc will suppose that n > m. The following result is due to Konig and 
O'Connell (2001), Katori and Tancmura (2004) and Bru in the real case (1989a) : 

Theorem 1. Let Ai(t),--- ,Xmit) denote the eigenvalues of Xf. Suppose that at time 
t = 0, all the eigenvalues are distinct. Then, the eigenvalue process {\i{t) , . . . , \m{t)) 
satisfies the following stochastic differential system: 



dXiit) = 2^/>^ dPiit) + 2 



^ Xiit) + Xk{t) 



dt 1 < i < m, t < T, 



where the {Pi)i<i<m o,re independent Brownian motions and t is defined by 
T := mi{t,Xi{t) = Xj{t) for some{i,j)}. 

Remark 5. With the help of the SDE satisfied by the eigenvalues, we can compute the 
ones satisfied by both processes (tr(X)) and (dct(X)): the former is done. For the 
second, we find that for t < Tq := inf{t,det{Xt) = 0} and for r G M: 

d{det{Xt)) = 2 det{Xt)^tT{Xf^)di^t + 2(n - m + 1) det{Xt) tr(Xj-^)dt 

d(log{det{Xt)) = 2^t^{X^)dut + 2(n - m) tT{Xf^)dt, 



d{det{XtY) = 2r{det{Xt)y ^ti{X^^)dut + 2r(n - m + r){det{Xt)Y tT:{X^^)dt 
so, we can see that for n = m, log(det(X) is a local martingale and so is (det(X))"^~". 

Lemma 1. Take Xq G Hm- Then for n > m, Xt G H^- 

Proof: In fact, this result is a direct consequence of the fact that for n = m, log det(X) 
is a local martingale, and so is (det(X))™~"'. Hence, for n > m, these two continuous 
processes tend to infinity when t ^ Tq which is possible only if Tq = oo, because every 
continuous local martingale is a time-changed Brownian motion. □ 

Corollary 1. // Ai(0) > • • • > Xm{0), then, the process U defined by 
is a local martingale. 

Proof : We could follow the proof given by Bru (1989a) or make straightforward 
computations using the derivatives of the Vandcrmonde function. But we prefer use a 
result from Konig and O'Connell (2001): for n > m, the eigenvalue process is the V- 
transform (in the Doob sense) of the process obtained from m independent BESQ{2{n — 
m + 1)). Thus, if G and G denote respectively the infinitesimal generators of these two 
processes, then, G{h) = and, for all function /, 

Gif) = l^G{Vf) ^ GiU) = 1g(1) = 0. 

Corollary 2. // at time t = 0, the eigenvalues of X are distinct, then, they will never 
collide, i.e. r = oo almost surely. 
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Proof : This result follows from the fact that the continuous process U tends to infinity 
when t ^ T which is possible only if r = cxd almost surely (We use the same argument 
as before). □ 

3.1. Additivity Property. The proof of this result is similar to the one derived by 
Bru in the real case (1989): 

Proposition 2. If{Xt)t>o oind iYt)t>o are two independent Laguerre processes L{n, m, Xq) 
and L{p, m, Yq) respectively, then the process {Xt + Yt)t>o is a Laguerre process L{n + 
p,m,Xo + Yq). 

Now, we introduce the Laguerre processes of noninteger dimensions 6. 

4. Laguerre Processes with noninteger dimensions 

Let X be a Laguerre process L{n, m, Xq) with n > m. If Xq G Hm, and if ^/Xt stands 
for the symmetric matrix square root of Xt , it is easy to show that the matrix O defined 
by Ot := ^/Tt'^Bf, where Xt = B^Bt, satisfies O^'O = OO* = Im- Thus, 



-"t - 

d-ft = OtdBt = y/xT^BtdBt 
is a m X m complex Brownian matrix. Replacing this expression in ([1]), one obtains : 

dXt = \fXtd-^t + dft^fXt^ Inl^dt 

Theorem 2. If {Bt) is a m x m complex Brownian matrix, then for every Xq G 
and for all 5 > m, the SDE 



dXt = VXtdBt + dBty'Xt + 25lrndt (4) 

has a unique strong solution in Hm- Furthermore, if the eigenvalues are distinct at time 
t = Q, then they satisfy the stochastic differential system: 



d\iit) = 2y%(t)d(3^{t) + 2 



^A,(t)-A,(t) 



dt 1 < i < m, 



where the {Pi)i<i<m are independent Brownian motions. 

Proof : The proof of the second part of the theorem is the same as before with S 
instead of n. So, we have to prove the first part. Note first that (det(Xt)), (logdet(Xt)) 
and {det{XtY) verify the same SDE with 5 instead of n. Hence, arguing as before, we 
can see that Tq = oo almost surely. On the other hand, the map a i— > a^^"^ is analytic 

in Hm (see Roger and Williams 1987, p. 134), so, the SDE has a unique strong solution 
for all t>0. □ 

Definition 2. Such a process is called the Laguerre process of dimension 6, size m and 
initial state Xq. It will be denoted by L{6,m, Xq). 

Remark 6. Any process {Xt)t>Q solution of (jl]) is a diffusion whose infinitesimal 
generator is given by: 

^ = 26 tr(K(^)) + 2[tr(xK(^)2) + tr(y9(^)2)] 



Remark 7. A simple computation shows that 

d{Xij,Xki)t = 2{Xii(t)dkj + Xkj{t)5u)dt, for all i,j,k,l G {1,... ,m} 

Now, we focus on both existence and uniqueness when 6 > m — 1 and Xq G (see 
Bru (1989) for the real case). 

4.1. The Process X~^ . If X is a Hermitian matrix, let X~^ be the Hermitian matrix 
max(X, 0). If we denote by (Aj) the eigenvalues of X, then (A^ = max(Ai,0)) are those 
of X+. 

Theorem 3. For all 5 £ and Xq = x £ Hm , the stochastic differential equation 

dXt = ^/x+dBt + dB^^fx+ + 26Imdt (5) 

has a solution in Hm- 

Proof : The mapping a \/a+ is continuous on Hm- Hence, X exists up to its 
explosion time (Ikeda and Watanabe 1989, theorem. 2. 3). Furthermore, from 

|2 I Ilrrll2^i:2 , \ \ v\\2 ^ /^/-i , ||x^||2\ 



||VX+||^ + <6^ + \\X\\^ < C(l + ||X| 

we can deduce that this explosion time is infinite almost surely (Ikeda and Watanabe 
1989, theorem 2. 4). ■ 

Proposition 3. // Ai(0) > ... > Am(0) > 0, then, for all t < S := inf{t, Aj = 
Xj for some (i, jj}, the eigenvalues of X^ verify the following differential system: 

dX.{t) = 2^fx[{t}di.,{t) + 2(5 + ^ X(f)-xl(t) I 1 ^ ^ ^ 

Proof: This differential system can be shown in the same way as in theorem [1] using : 

{dX,j,dXki)t = 2{Xf^{t)5kj + X+{t)dii)dt, for alH, j, k,l £ {I, . . . ,m} 

Proposition 4. // Ai(0) > . . . > Am(0) > 0, then, for all 6 > m - 1, t > 0, Xm{t) > . 

Proof : First, we note that S = 00 almost surely. Indeed, one can easily show that 
the process U defined by : 

U{Xi{t),...,Xm{t)) ^ 



n.<,(AiW-Ai(t)) 

is a local martingale. For the proof, we follow in the same way as Bru (1989). 

Theorem 4. // Ai(0) > ... > Am,(0) > 0, then, for all 6 > m — 1, Q has a unique 
solution in in the sense of probability law. 

Proof : By Proposition HI the solution of the SDE (jS]) remains positive for all t > 0, 
thus, it is a solution of ([5]). □ 

Theorem 5. Let be the space of positive Hermitian matrices. Then, whenever 
the SDE ^ has a solution in H^, for fixed t, its distribution is given by its Laplace 
transform: 

Exo(exp-(trnXt)) = (det(/™ + ^to))-'^ exp(-tr(Xo(/„ + 2tu)-^u)), (6) 
for all u in 



Proof : For s G let g{t, s) = A^^ exp(-F(t, s)) , where 

At = det(/„ + 2ut), Wt = {Im + 2ut)-^ u, V{t, s) = tv{sWt), 

First, note that W € Hm- To proceed, we need a lemma : 

dg 

Lemma 2. The function g satisfies the heat equation: — = ^g where =Sf is the infin- 
itesimal generator of X. 

Proof of the lemma : If we write s = x + iy, then, using the fact that x is symmetric, y 
is skew-symmetric and W is Hermitian, we can see that tr(sWt) = tr(a;M + iyN) where 

W + W W -W 

M = N = - —. 

2 2 

Observing that = M and A''^ = —A'', we can deduce that g satisfies the conditions 
of Proposition 1. Besides, 

^ = -g{25iT{Wt)-2iT{sW^)) 



dxdy dydx , 
tr(x(^0-0)) = 5tr(x(M2+Ar2))=5tr(xW2). 

Finally, noting that tr(M) = tr(W^), we obtain the equality. Now, we consider the 
process {Z{t^ -^t)) defined by X^) = g{t\ — t, Xt) for all t < ti for fixed ti. From the 
lemma, we deduce that Z is a bounded local martingale and thus is a martingale. So, 
the result follows from a simple application of the optional stopping theorem. ■ 

Corollary 3. Let {Xt)t>o be a Laguerre process L{S,m,x) where x G Hm^ ■ For 6 > 
m — 1, its semi-group is given by the following density: 

with respect to Lebesgue measure dy = Y[p<q dy^g Y[p<q dypg where y = y^ + iy'^ and qFi 
is a hypergeometric function of Hermitian matrix argument (Chikuze 1976, Gross and 
Richards 1991). 

Proof : In fact, this result can be easily deduced from the case where 5 = n\s integer, 
since, in this case, Xt is a non-central complex Wishart variable W{n,2tlm,x) (James 
1964) with density given by: 

^ (2t)^"r^(n) ^""P ^'^"^ + ^^^^ yT~""oFiin;^)l{y:^o} 



with respect to dy. Hence, taking 6 instead of n and denoting by Wt this new variable 
(starting from x), we can see that : (we wih use \y\ to denote det(y)) 



^.(e— ) - (^^yjf^e-^ / exp (-ltr((/ + 2n.).)) |.|-V.(.;g 



)dy 



ly>0 

2^»Ti(5|^|-<5 t f 1 1 

e"^ / exp{-2ttv{x^2(I + 2ut)x^2z)\z\^-"'oFi{5;z)dz 

Jz>0 



= exp(-— )|/ + 2ut\~' exp (tr(^(/ + 2ut)-') 

= |/ + 2tit|"^exp ^-^tr(x(/ + 2ut)"^(/ + 2iit-/))^ 

= |/ + 2ut\~^ exp (- tr(x(/ + 2ut)~^u)) 
which is equal to (l6|). ■ 

Remark 8. In the last proof, we used the change of variables z = x^/^yx^/^ which 
gives dz = \x\"''dy. For the second integral, see Faraut and Koranyi (1994), proposition 
XV.1.3. 

Remark 9. The expression of the semi-group extends continuously to the degenerate 
case, namely: 

where 0^ denotes the null matrix. 

Corollary 4. For 6 > m — 1, the semi-group of eigenvalue process is given by: 



V{x) \2t \x 

where x = (xi, . . .,Xm),y = {yi, ■■■,ym) so that xi > . . . > Xm > 0, yi > . . . > ym > 0, 
5 = m + v such that u > —1 and ly denotes the modified Bess el function (Lebedev 1972). 

Proof : The expression of the semi-group can be computed using Karlin and MacGre- 
gor formula (1959) since, for 6 > m — 1, the eigenvalue process is the /i-transform of the 
process consisted of m independent BESQ{2(5 — m-\-l)) conditioned never to collide, as 
stated by Konig and O'Connell (2001). Another proof is given by Peche (2003, p. 68). 
Here, we will deduce the expression of qt{x,y) from pt{x,y) following Muirhead (1982), 
namely, by projection on the unitary group : we will use Weyl integration formula, then 
give a determinantal representation of hypergeometric functions of two matrix argu- 
ments. First, we state Weyl integration formula (Faraut 2006) in the complex case: for 
any Borel function /, 

/ f{A)dA = Cm / f{uau*)a{du){V{a))'^dai . . . dam, 

JHm Ju{m) Jr™ 

^m(m-l) 

where Cm = — ^, U{m) is the unitary group, a is the normalized Haar measure on 

Tm{m) 

U{m), a = diag{ai) and A = uau* . Hence, the semi-group of the eigenvalue process is 



given by (James 1964): 

qtix,y) = Cm{V(yf) / pt{x,uyu*)a{du) 
Ju{m) 

Cm{V{yf) fV _^i±M. f-Pf V /" TT ix ^uyu* 



n'"'"'-"{V[,,f) -ft fill! /ft i -^ 

= (2t)".(".«)r„(™)r„(m + .)n; " (n»j o^-iC^ + ^ij.;!/). 

where y = diagiyj) , x is a positive definite matrix whose eigenvalues are xi, . . . ,0;^, 
o-Fi in the third hne is an hypergeometric function with two matrix arguments (Gross 
and Richards, 1991) and 5 = m + v,i' > —1. Next, we need a lemma. 

Lemma 3. Let B,C E and let (6^), (cj) be respectively their eigenvalues. Then, 

^{m-i), -,^ -A- frfu,- + 1))"* 



n 



Tmim + det (p^qiim + l)i<i<p, (1 + (pj)i<j<q; hcf)ij 



il(r(<^, + l))- V{B)V{C) 

for all > —1, l<i<p, l<j<q- 

Proof : Recall that the hypergeometric function of two matrix arguments is given by 
the following series: 

F((n^ -p.r^ vv nr=iK)ra(i?)a(c) 

pi^,((a,)i<,<p,(e,)i<,-<„B,C) _ 2_^2. f^^^^— c;(7)^' 



It is well known that 



k\dr 

Cr{B) = — — ^Sr(6l,...,6m), 

[m)r 



where is the Schur function and di- = 5^(1, • • • , 1) is the representation trace or degree 
(Gross and Richards 1991, Faraut 2006). Substituting in the series, one gets: 

Mim + M,)i<Kp, {rn + </,,)i<,<„ B,C)-^l^ \{'=M + <^,)r (m). ' 

/c — ^ 

Now, we write: 

. _ -^1- r(^i + m + kr -r + l) r{fii + 1 + kr + 5r) 

[m + fii)r-ll r(^, + m-r + l) " ii r(/.i + m - r + 1) 
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Doing the same for each (m + 4'j)T and for {m)-^, we can see that: 

? — 1 .7 — 1 



fc=0 r r= 



/ ?n — 1) \ 

where /3 = ( 2 ) iP~i~^)- S^* desired result, we use the "Hua formula" 

(Faraut 2006) : 

Lemma 4. Given an entire function f, i.e. f{z) = ^1^=0^^^^ , then 



Thus, we get: 



V{B)V{C) y s^ilrr 



n 



det {p^q{{fJ,i + l)l<j<p, (1 + (t>j)l<j<q; blCf)if 

V{B)V{C) 

For p = and g > 1, it reads : 

oF,((m + 0,).<,<,; B. C) = .-=^(«+'>r„(™) n 

det (o^q((l + <^j)i<j<q; 
y(B)F(C) 

and similarly, 

° °^ ^^^^ 

which can be viewed as Harish-Chandra formula for the "Itzykson-Zuber" integral (Collins 

X 

2003). We now proceed to the end of the proof. TaJiing p = 0, q = 1, B = C = y, 
we get: 

p . , A ~^ = (4t2)Mrn-i)/2r^(^ + ^)r^(^) det(o^i((i/ + l);^)) 
oi^i[m + u, ^^^,y) W"^-i)(r(zy+l))"^ V{x)V{y) 
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The expression of qt{x,y) follows from a simple computation and from the fact that: 
o^l{{u + l)■,x^y,/4e))iJ _ f 2t y 

r(i. + i) \vmJ t 

Proposition 5. The measure defined by p{dx) = (det(x))''~'" dx on Hm^ is invariant 
under the semi-group, i. e, pPt = p. 

Proof. Denote by Pt the semi-group of Laguerre process L{6,m,x) for 5 > m — 1. 
Then, we have to show that 

/ Ptf{x)p{dx) = [ f{y)p{dy) f e Co{hJ). 

Jx>0 Jy>0 

This follows by a similar computation and the same arguments as in the proof of corollary 

m ■ 

Remark 10. For Wishart processes, it is easy to see that p{dx) := (det(x)) 2 ^ l^^yQ^^dx 
is invariant under the semi-group. 

5. GiRSANOv Formula and Absolute-continuity Relations 

The index v > —1 of a L{6,m,x) is defined hy 1/ = 6 — m. In this section, we will 
discuss in the same way as in [] to derive absolute-continuity relations between different 
indices. 

5.1. Positive Indices. Take a matrix-valued Hermitian predictable process H. Let 
be the probability law of L{6, m, x) for 6 > m — 1 and x G Hm- Define: 

'•^ ti{H,dB, + TTsdB's) 



1 r 

$i = exp(Lt-- / tr(i/2)ds), 

^ JO 



where i? is a complex Brownian matrix under Q^. We can easily see that the process (5 
defined by (3t = Bt — Jq Hgds is a Brownian matrix under the probability 

Pfl^, := ^t-QiW,, 

Furthermore, {Xt)t>o is a solution of 

dXt = y%dl3t + d(5ty/Xt + (VYtHt + Hty/X't + 25Im)dt. (7) 

For Ht = fy/Xl ^, ([7]) becomes 

dXt = y/Xtdpt + dptVx't + 2{6 + i^)Irndt, 

so that {Xt)t>o is a L{6 -\- iy,m,x) under . Thus, we proved that : 

Theorem 6. For 6 > m — 1, 

= exp j\T{y/xr'dBs + ^/x's'^dB,) - ^ j\T{Xs-^)d^ ■ QIW,. (8) 
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Proposition 6. 



(9) 



Proof : We know that V„(det(n)) = det(n)u~^, hence, V„(log(det(u))) = Then, 
using the fact that for 6 = m, (log(det(Xt))) is a local martingale, we get from Ito 
formula that: 

log(det(Xt)) = log(det(Xo)) + f ti{X-\^sdBs + dB*^s)) 

Jo 

= log(det(Xo)) + / tl{^/x's'^dBs + ^/Ys'^dBs). □ 
Jo 

From ([9|), it follows that: 
Corollary 5. 

f\ ^^-U. A ,^ ^ det(y)-^/2p™+'^(x,y) 



Q™(exp {^-yI tr(^7')rfs) \Xt = y) 



det(x) pY^{x,y) 



-^f^(det(z))^/^°M!l±i^ 



where z = — -jr . 

Now, we state the following asymptotic result: 
Corollary 6. Let X be a Laguerre process L(m,m,x), then, as t — > cxd; 

{mlogty Jq 

where Tj is the first hitting time of 1 by a standard Brownian motion (3. 
Proof : From ([9]), we deduce that: 

oFi{m + 2z^/mlogt,xy/4t^) 
oFi(m,xy/4t2) 

Noting that (^^™-)-'^/»"log* = g-f ^ ^^d since both hyper geometric functions converge to 1 
as t — > oo, we obtain: 



t~*oo 



e 



-V 
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Then, since 



lim t^'pTix,2y) = lim £_!^e-*'^(^)oFi(m, ^) 

p- tr(y) 



r^(m) 

we get: 



= f Q7 (exp(- . , / tr(JC7l)ds|X, = (») f'pT{x,ty)dy 

Jy>0 \ {mlOgl) Jq J 

— ^ e , 

by dominated convergence Theorem. ■ 

5.2. Negative Indices. Take < a < det(a;). The same computation as in parag. 5.1 
with Hf = —vy/Xt ^ ,0 < v < I, shows that 

where Ta := inf{t, det(Xt) = a}. Letting a — > and using the fact that Tq = oo a.s 
under Q™, we get : 

_ f det{x) Y ^rn+uj 

Proposition 7. For allt > and < < 1, 

Proof : Prom the absolute-continuity relation above, wc deduce that : 
On the other hand, using the expression of the semi-group, one has : 

by Kummer relation (cf Th 7. 4. 3 in Muirhead 1982). Taking s = —u, we are done. ■ 

13 



6. Generalized Hartman- Watson law 

Henceforth, we will write ^ to denote one-dimensional hyper geometric functions. We 
define the generalized Hartman- Watson law as the law of 



[ tv{X-^)ds under Q'^{-\Xt = y). 
Jo 



Its Laplace transform is given by: 

Qx exp — - / tr(X, )ds \Xt = y) = — — -det(z / — — (10 

\ 2 Jq J Tm{m + u) oFi{m,z) 

z = xy/4:t'^. Recall that for m = 1, this is the well-known Hartman- Watson law and 
that its density was computed by Yor (1980). Here, we will investigate the case m = 2. 
The Gross and Richards formula is written for p = and q = 1 : 

j:. , , N dQi{z^~\^i{m + v -j + 1, Zi)) 
oF,{m + u,z) = , 

where (zj) denote the eigenvalues of z and V{z) = Yli^j{zi — Zj) is the Vandermonde 
determinant . Noting that Tmi'rn + v) = YYjLi -^("^ + ~ J + 1); then : 

det(zf-^')/^W,-,(2^)) 

^ det{zf^/'lm-j{2^)) 
Without loss of generality, we will take t = 1. 

Proposition 8. For m = 2, let Ai > A2 be the eigenvalues of ^/xy. Then, the density 
of the generalized Hartman- Watson law is given by : 

^_ ^/hMv zsinh(p^/r^^)e-^v'^^^°'^^^^e- ''"^""'' (sinh|/) sin{^)dzdy 



pvr^/W /o /o "'°'''^^!^3^/o(gVA7A^^x)rfnrfx 

for u > 0, where p = Ai — A2. Furthermore, if Ai = A2 := A, then: 

AXve^ g{y)e~^{suih.y)sm.{^)dy 



^^"^ 7r2^/W i^2(^;l;2;A2) 
where 

1 vr /2(2A cosh y) -|- L2(2A cosh y) 
9[y) = 3 + 2 2Acoshy ' 

and L2 is the Struve function (Gradshteyn and Ryzhik, 1994)- 

Proof : For m = 2 , (llOp becomes: 

^ Ai4+i(Ai)/.(A2) - A2J.+i(A2)/.(Ai) 

^ Ai/;(A,)/o(A2)-A2/,(A2)/o(A;) ' 

so, using the integral representations below (Brychkov, Marichev, Prudnikov 1986, p. 
46) : 



x{aljyj^i{ax)ly{bx) — bljjj^i{bx)ly{ax)) = {a —b)j uly{au)ly{bu)du 
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X 







with X = l,a = Xi,b = A2, and (Gradshteyn and Ryzhik 1994, p. 734): 



where a > 0, 3ft(z/) > —1, with a = 1,6 = (Ai — \2)u := pu et c = 2\/ A1A2U, the 
numerator of (fTO]l is then equal to: 



Taking 1/ = 0, the denominator is then equal to: 

^ Jo Jo V 1 - a;"' 

Thus, pO]) becomes: 

£ /o "-MWT^) /,,(gVA7A^^x)^^^:^ 



/o /o --';(^^) /,(gVA7A^^x)rfz.rfx 



Now, we only have to use the integral representation of (Yor 1980): 

jj2\/ Xi\2ux cosh u> —2uuj 



2z7r Jc 



c 



2ivr 7c; 7o (27ru3)i/2 

where C is the contour indicated in Yor (1980), hence, the density function is given by: 



^ 1 /o /o w -'-^^g^^ e^v^— ^-e-^^ndxd^ _^ 
ivr^/W /o Jo "'°^^^';"^I^/Q(gVA7A^nx)rft^rfx ^'^"^ 

We can simplify this expression by integrating over C to see that the numerator is equal 
to (Yor 1980): 

ttV2itv^ Jo Jo Jo VI -x2 V 

Setting z = ux, The numerator is written 



ttV2ttv^ Jo Jo Jo yu^ — z"^ 

that we can integrate with respect to u to get 

p7rV2TTv^ Jo Jo 

Now, we prove the second part. In this case, p = and we have to evaluate : 
A^e^ Jo7;/o°°^e-^^— t-^e-^(sinh^)sin(^)rfndxdy 

loJo^'^^dudx 
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Setting z = ux, the numerator reads : 




^0 



Integration with respect to z yields (Gradshteyn and Ryzhik 1994, p. 369): 

^ no 

Xve V f . Attv 

g{y)e - {shy)sm{ )dy 



'K\p2/JVlfi Jo V 

For the denominator, we use the fact that ■j^{zli{z)) = zIo{z), which yields: 

1 '■^ uIo(2Xux) , , Ii(2Xx) , 
duax = / „ ax 




/o Vl — Jo 2Ax\/T — x^ 

Then, the following formula 

' x'^-^ic? - x'^)^-'^I (cx)dx - 2-^-ia2/^+°+'^-2c^ r(/?)r((a + ;y)/2) 

X [a X) I,[cx)dx-Z a T{P + {a + u)/2)T{u + 1) 

^ ,a + V ^ a + u ^ cP'c? , 

i^2(^-; /3 + + 1; ^) 

taken with a = 0, a = 1, /3 = 1/2, c = 2A, v = \ gives: 



We can proceed differently : let Ai = A2 + /i then pUj) reads: 

((A2 + ; (A2 + /l)/.(A2) - A2/.+i (A2)4(A2 + 

((A2 + K)h (A2 + /i)/o(A2) - A2/i (A2)/o(A2 + h))lh ■ 

Next, we let /i — > 0. As usual, we first compute the numerator and then take 1^ = 0. To 
do this, we shall evaluate : 

^ ^ lini + ^^^'^+^ + ~ ^^^'^+^ ^^^^ 

4(A2 + M - /.(A2) 
ij = lim 

h^O h 

which are equal respectively to -^{xl^^i {x)) and taken for x = A = Ai = A2. 

dx ^ dx 

Using the differentiation formula — (x'^/j^fx)) = x^Iy^iix) (Lebedev 1972, p. 110), we 

dx 

get: 

-^{xlu+i (x)) = -vlv+i {x) + x/;,(x), -^{Iu{x)) = --Iu{x) + ly^l (x), 
dx dx X 



thus: 



N = /,(A)(-i./,+;(A) + A4(A))-A/,+;(A)(-^4(A) + 4_,(A)) 



A 



A(42(A)-4+,(A)4_,(A)) 
4^(A)-4+j(A)4_j(A) 
/o2(A)-/;(A)/_,(A) 
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Using the integral representation below (Gradshteyn and Ryzhik 1994, p. 757): 

I^,{z)U{z) = -[ ' cosiifi - 1^)6) I^+^{2z cos 9)d9, + i/) > -1. 

^ Jo 

the numerator is written as : 

2 r^/^ 4 /-Va 



N 



-f {I -cos 29) I2^{2X cos 9)de = - f {sin^ 9)1 2u{2\ cos 9)de 
Jo ^ Jo 

= - [ Vl-r'^l2u{2Xr)dr, 
71" Jo 

Thus, using (|lip . the denominator is equal to 

D = ^ /'Vl-r2/o(2Ar)dr = ^i^2(J;1;2;A2) 

Jo 4 2 

Finally, the integral representation of gives : 

Xue'^^Vn /o°° sinh(y) sin (^^^ r^T^e-^^'^^'^i^^'dr citx 



7rV27ru3 /gVl - r'^Io{2Xr)dr 



TT 



V27r7x3 /g^^r^72/o(2Ar)dr 
7. The Law of Tq 

Recall that: For < < 1, 

Qr-in > t) = r (m) det(^).^Fi(^,m + ^,-^) 

+ 2r 2r 

Proposition 9. Lei m = 2 and Xi > X2 be the eigenvalues of x. The density of 
So := 1/{2Tq) under Q™"*^ is given by: 

. ^ (AiA2)"^^"-'e^(^i+^^)» 1^1(2, u + Am) - 1^1(2, u + l, X2U) 

r(z. + i)r(z.) (A1-A2) 

Corollary 7. // Ai = A2 := A, the density is written: 

V[y + 2)\:[y) 

Proof : Recall first that when m = 1, So = Ju/x, where 7j/ is a Gamma variable 
with density r'^~^e~^dr. With the help of the Gross-Richards formula, it follows that 
for m = 2, 

QT"{So <u) = T^T^^T— T^^^'"(Aii^i(z^, 1^ + 2, -Ain)i^i(z. - 1, + 1, -Aau) 

- A21=^l(l^, U + 2, -X2u)i^i{v - 1, 1/ + 1, -Alii)), 
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where 5*0 := l/(2To). This is a C°° function in u. Hence, we wih compute its derivative 
to get the density. Recall that : 

^ 1^1 (a, b,z) = ^ 1=^1 (a + 1, 6 + 1, z), 

thus : 

fiu) = ^QT^iSo <u) = K{u,Xi,\2)u'"'-\A-B) 



where 

Kiv, Ai, A2 



(A1A2)' 



r2(i/ + 2)(Ai - A2) 

A = 2u{{Xii^i{u,iy + 2,-Xiu)i^i{i^-l,u + l,-X2u) 

- A2i^i(j^, iy + 2, -X2u)i^i{v - 1, 1/ + 1, -Xiu)) 

B = -^((A?ui^i(i/ + l,z/ + 3,-Aiu)i=^i(i/-l,z/ + l,-A2u) 

- Xlui^i{v + 1, 1/ + 3, -X2u)i^i{v - 1, 1/ + 1, -Am)). 
Then, we use the contiguous relation : 

hi^i{a,h,z) - 6 1=^1 (a - l,6,z) = z i=^i(a, 6 + 1, z) 

to see that 

Aini^i(z^+ + 3,-Aiu) = {u + 2){i^i{v,v + 2,-Xiu) - i^i{v + l,v + 2,-Xiu)) 

X2Ui^i{u + + ^,-X2u) = (i/ + 2)(i#"i(i/,z/ + 2,-A2m) - + l,v + 2,-X2u)) 

implies that: 

f{u) = Ki{u, Xi,X2W-HC + D-E-F) 

where 

K(.X - KA1A2)- 

- r2(. + 2)(A,-A2) 
C = Aii^i(z/, 1/ + 2, -Ain)i^i(z/- l,zy + 1,-A2n) 

= Aii,Fi(zy + l,z^ + 2,-Ai77,)i,Fi(zy- l,zy + l,-A2«) 
E = A2i^i(/^,i^ + 2,-A2n)iJ^i(zy- l,z^ + l,-Ain) 
F = A2i^i(/> + l,J^ + 2, -A2n)i^i(z/- + 1,-Ai«), 

Applying again the above contiguous relation yields: 
XiUi^i{u + l,iy + 2,-Xiu) = + 1) (1^1(1/, Z/ + 1,- Am) - i^i(z/ + + 1,-Am)) 
X2Ul^l{l' + l,u + 2, -X2U) = {u + -X2U) - i^i{u + 1, 1/ + 1, -A2^i)) 

X2Ui^i{u, u + 2, -X2U) = {u + - 1, 1/ + 1, -A2tx) - 1=^1(1^, + -X2U)) 

Xiui^i{u, u + 2, -Am) = {u + - 1, 1/ + 1, -Am) - v + l, -Am)) 

Replacing in the expression of /, we obtain 

f{u) = K2{v,XiA2)u''''-\G - H), 
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where 



K2{v,\l,\2) = 



T2{u + 2)i\i- \2) 
G = + + l,-\2u)i^i{v + l,-Xiu) 

H = i^i{u + l,u + l,-Xiu)i^i{u-l,u+l,-X2u) 

Eventually, writing 

r2{iy + 2) = r(z^ + 2)r(z/ + 1) = 1/(1/ + i)r(i/ + i)r(z/) 

i^i{a,a,z) = i^i{a,b,—z)=e'~^i^i{b — a,b,z), 

we get 

(AiA2)^u2;.-2e-(^i+^2)n ^^^(2, + 1, Aiw) - i^i(2, 1/ + 1, Agu) 



fiu) 



r(i/+i)rM (A1-A2) 

Th case Ai = A2 is treated in the same way as before (for the Hartman- Watson law). In 
fact , writing Ai = A2 + ^ and letting h ^ 0, we see that the density is given by : 

2\2v y2u-l ^-\u 

= r(. + 2)r(.) 1^1(^-1^- + 2, -A.). 

8. Conclusion 

The Gross-Richards formula has been the main ingredient in this paper, since it 
enables us to express more explicitly the special functions of matrix argument. The case 
m = 3 can be treated in the same way, but computation becomes too complicated. So, 
if we want to deal with the general case, it will be convenient to find a more explicit 
formula. Indeed, Schur functions can be expressed as polynomials in the elementary 
symmetric functions or as polynomials in the completely symmetric functions hr- 
More precisely, we have : 

sx = det(eA,_i+j) 1 < ij < n 
sx = det{hx'^_i+j) l<i,j<n 

where A is a partition of length < n, and A' is the conjugate of A (Macdonald 1995). So, 
using these two identities, can we improve our results? 

9. Appendix: special functions 

9.1. The hypergeometric series. The multivariate hyper geometric functions were 
studied by Muirhead (1982) in the real symmetric case, Chikuze (1976) for the complex 
Hermitian case and Faraut and Koranyi (1994) in a more general setting. For Hermitian 
matrix argument, they are defined by: 

pFq{{ai)i<i<p, {bj)i<j<q; X) = J2Y1 TaT ThY~~T^ 
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where r = (ki, . . . , km) is a partition of weight k and length m such that ki > . . . > km, 
(o)t- is the generahsed Pochammer symbol defined by: 

, . ^ T{a + ki-i + l) 

i^)r = ll r(a-Z + l) ' r={k„...,km) 

1=1 

and Ct-{X) is the zonal polynomial of X such that : 

T±k 

Several normalizations for this polynomial exist in the litteraturc but we consider this 
one. This polynomial is symmetric, homogeneous, of degree k in the eigenvalues of X 
and is an eigenfunction of the following differential operator : 



1 = 1 ' 1 = 1 l<feytt<m 

Besides, it is identified with the Schur function Sr and Cr{YX) = Cr{VYXVY) for any 
Hermitian matrix Y. It is well-known that, p = q + 1, then the hyper geometric series 
is convergent for < ||X|| < 1 (j| • || is the norm given by the spectral radius) , ii p < q, 
then it converges everywhere and else, it diverges. 

9.2. The modified Bessel function (Lebedev 1972). The modified Bessel function 
with index G R is given by the following series : 

'-i^^^E kWl. + k + D (2) ' 

It can be represented through standard hypergeometric functions 0^1 and iFi : 
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